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ABSTRACT 



We use ba r yon ac oustic peak measurements by lEisenstein et al.l (120051 ) and 



Percival et al.l (j2007al ). together with the WMAP measurement of the apparent 



acoustic ho r izon a ngle, and galaxy cluster gas mass fraction measurements of 



Allen et al.l (120081 ). to constrain a slowly rolling scalar field dark energy model, 



0CDM, in which dark energy's energy density changes in time. We also compare 
our 0CDM results with those derived for two more common dark energy mod- 
els: the time-independent cosmological constant model, ACDM, and the XCDM 
parametriz ation of dark energy's equation of state. For time-independent dark 
energy, the lPercival et al.l (j2007al ) measurements effectively constrain spatial cur- 
vature and favor a close to spatially-flat model, mostly d ue to the WMAP CMB 
prior used in the analysis. In a spatially-flat model the iPercival et al.l (j2007al ) 
data less effectively constrain time-varying dark energy. The joint baryon acous- 
tic peak and galaxy cluster gas mass constraints on 0CDM model are consistent 
with but tighter than those derived from other data. A time- independent cos- 
mological constant in a spatially-flat model provides a good fit to the joint data, 
while the a parameter in the inverse power law potential 0CDM model is con- 
strained to be less than about 4 at 3a confidence level. 
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1. Introduction 

About a decade ago type la supernova (S NIa) observations provided initial evide nce 
that the cosmological expansion is accelerating (IRiess et al.lll998l ; iPerlmutter et al.lll999l ). If 
general relativity is valid on the scales of current cosmological observations, more recent SNIa 
data as well as results of various other cosmological tests including large-scale structure tests 
and observations of the cosmic microwave background (CMB) anisotropy can be reasonably 
well reconciled if we assume that about two-thirds of the cosmological energy budget is in 
the form of dark energy. 

Different theoretical models of dark energy have been proposed over the years. A few 
models try to do away with the need for an exotic dark energy component by modifying gen- 



eral relativity on large scales (see, e.g.. IWang 



2008 



Demianski et al.ll2008l ; iTsujikawa et al. 



20081 : ICapozziello et al.l bood : fTOl20od : lGannouii fe Polarskil booi T If general relativity 



is 



valid we need a substance that has negative pressure, p < —p/3 (where p is the energy 
density), to have accelerated c osmological expansion. The simplest standard cosmological 
model is ACDM ([Peebles! 1 19841 ) in which the cosmological constant A has negative pressure 
and powers the current accelerated expansion of the universe. Although A has a quantum 
field theory motivation as vacuum energy, ACDM has a number of apparent problems. The 
most celebrated is the fact that the value of vacuum energy density calculated from field the- 
ory with a Planck scale cutoff is many orders of magnitude larger than the measured value. 
Because of this other models have been developed, despite the fact that the simple ACDM 
model provides a fairly good fit to most cosmologic al data. In our paper we also study the 
slowl y rolling scalar field dark energy model (0CDM, iPeebles fc Ratralll988l ; iRatra fc Peebles 
19881 ). In the 0CDM model the small (classical) value of the current vacuum energy density 
is a consequence of the scalar field dynamics. The third model we consider is the XCDM 
parametrization. XCDM parametrizes dark energy's equation of state as p x = u x p x , where 
uo x is a negative co nstant. This approximation is not accurate in the scalar field dominated 
epoch (IRatra SEE 



In the 0CDM model one can explain the accelerated expansion of the universe by in- 
troducing a scalar field minimally coupled to gravity. The action for such a term is 



1 For recent rev iews of dark energy see, e.g., Ratra fe Vogelev ( 20081 ). Linder ( 2008 ). Frieman et al. ( 2008 ). 



and Martin 



(200 



We assume that dark ener gy and dark matt er only couple gravitationally. For discussion of model s 



with other couplin gs see, e.g., 



Costa et al 



He fe Wand (hoosl). and blivares et alj (boosl). For oth e r models of d ark 



( 20081 ). iMainini fc Bonomettol (j2007t ). iBrookfield et all (12 008). 



energy see, e.g. 



Grande ct al. 

(|2007t ). lNeupane fc Schererl (|2008l ). lMathews et all (|2008f l. lUsmani et al.l (|2008h . and llchiki fc Keuml |2008L 
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S,k 



d 4 Xy/—g 



G 



where G is the gravitational constant. If the scalar field is close to homogeneous on cosmo- 
logical scales, then to leading order it's energy density and pressure are given by 



1 fd^y v{4>) 

2 left / G 



(2) 



dt 



V(4>) 
G 



(3) 



When the scalar field changes only slowly in time, the effective equation of state parameter 
ujcj, =V4>I P4> is negative and the scalar field acts like a time-dependent cosmological constant. 
To specify the 0CDM model one has to pick a specific form of potential energy density 
V ((/)). Neither cosmological observations nor fundamental particle physics theory can provide 
significant motivation for a specific form of potential energy and a lot of different cases 
have been studied. In our paper we work with the inverse power law potential energy 
density V(<p) oc 0~ a , because it has been well studied and it provides a practical way of 
parametrizing the slowly rolling scalar field with one nonnegative dimensionless parameter 
a. Physically, large values of a correspond to ra pid time evolution , while the limit of a = 
gives a time- independent cosmological constant. iPodariu &: Ratral (120001 . Fig. 2) relate this 
</>CDM model to the XCDM parametrization and discuss how a and effective are related. 
For large values of a the time-dependent equation of state parameter changes very fast 
and the XCDM parametrization fails to provide a good phenomenological description of the 
scalar field. Figure 1 shows the residuals between comoving distance calculated in 0CDM 
and XCDM models. Already for a = 2.0 the predictions of XCDM differ significantly at 
high redshifts. If a is very close to zero the scalar field equation of state changes very 
slowly and becomes more and more difficult to distinguish from ACDM and the XCDM 
parametrization is reasonable. If a turns out to be a very small but nonzero number, a 
lot of independent high precision cosmological measurements supplemented with the better 
understanding of underlying high energy physics will be necessary to discriminate between 
different dark energy models. 

Assuming the cold dark matter (CDM ) model of structure fo rmation (for a discussion 
of apparent problems with this model see iPeebles fc Ratral 120031 . and references therein) , 
and assuming that the dark energy is a time-independent cosmological constant (see, e.g ., 
Wang fc Mukherieel hooi bong et aJliooi llchikawa fc Takahashil |20o3 Ivirev et allbooa ). 
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CM B anisotropy data com bined with independent dark matter density measure ments (see, 



e.g.JChen fc Ratra 



2001bl ; IPage et al. 



2003 



2003b ) are consistent with negligible spatial c urvature (see, e.g. JPodariu et al. 



Spergel et al.l 120071 ; iDoran et al.l 120071 ). CMB anisotropy data in 



combination with the low measured density of nonrelativistic matter then require the pres- 
ence of dark energy and so are consistent with the SNIa results. 

Many different observational tests have been used to constrain cosmological parameters. 
An issue of great current interest is whether dark energy is Einstein's cosmological constant 
or whether it evolves slowly i n time and varies weakly in space. Current SNIa data are 
unable to resolve this (see, e.g., Mignone &: Bartelmann 2008 ; Wu fc Yull2008 ; Lin et al. 2008 



Dev et al.ll2008l ; iLiu et al.ll2008l; iKowalski et al.ll2008f ). but future SNIa data will improve the 
constraints (see, e.g., iPodariu et al.l l2001aj ) and, unless a has a very small value, might be 
able to detect time variation of dark energy. Current SNIa and CMB data are consistent 
with the ACDM model, but it is not yet possible to r eject other dark energy models with h igh 



statistical confidence (see, e.g.. iRapetti et al.ll2005l ; IWilson et al.ll2006l ; iDavis et al.l 120071 ) . 



To tighten the constraints on cosmological parameters, it is important to have many 
independent tests of dark energy models. Comparison of constraints from different tests 
can help uncover unknown systematic effects, and combinations of constraints from differ- 
ent tests can better discriminate between models. Other observational tests under recent 
discussio n include the angular size of radio source s and quasars as a function of redshif t 



Chen fc Rat~rall2003al :I Podariu et al1l2003l: balv et alJl2007t ISantos fc Lima 



Zhu & Sereno 


20 


osf), \ 


2008 




Dore et al. 


2007 



,see, e.g. 

strong gravitational lensing (see, e.g.. iLee &; Ng 



2007 



Qguri et al.ll2008l: IZhang et al 



0, 



2008 



2007; 



2008h. weak gravitational lensing (se e, e.g., iTakada & Bridle! 120071 : iFu et al.l 



2007; La Vacca k Co 



ter as a function of redshift (see, e.g. . 



ombol 120081). measurements of the Hubble parame 



Samushia fc Ratral 120061 ; lLazkoz fc Majerottd 12007 



Wei fc Zhang! 120081 ; ISzydlowski et al.ll2008r). large- s cale structure baryon acoustic oscilla- 
tion peak measurements (see, e.g., Kia et al.l 120071 ; iLima et al.l 120071 : ISapone fc Amendola 



2007 ) and galaxy cluster gas m ass fraction versus redshift data (see, e.g. 



Allen et al. 



2004 



Chen fc Ratral 12004!: ISenl 120081) . For recent reviews of the observational constraints on dark 



energy see, e.g. 



Kurek fc Szvdlowskvi (120081 ) and IWangl ((20071) 



Many different observational test have been used to constrain the slowly-rolling scalar 
field dark energy model. The constraints are getting tighter as the quality and quantity of 
new measurements is increasing. Constraints on the a parameter from different tests are 
shown in Table 1. In our paper we use baryon acoustic oscillations (BAO) peak measurements 
to constrain the 0CDM model of dark energy and compare our results to the constraints 
on ACDM model and XCDM paramet rization. Since the peak has been measured at only 
two redshifts, z = 0.2 and z = 0.35 (jEisenstein et al.l 120051 ; iPercival et al.l l2007al ). BAO 
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data alone can not tightly constrain the models. To more tightly constrain the dark energy 
models, we perform a joint a nalysis of the BA O data with new galaxy cluster gas mass 
fraction versus redshift data ( Allen et al. 20081 ) PI The resulting constraints are consistent 
with, but typically more constraining than, those derived from other data (see Table 1). 

In Sec. 2 we briefly describe the BAO method we use. In Sec. 3 we summarize the BAO 
and galaxy cluster gas mass fraction data and computations. We discuss our results in Sec. 
4. 



2. Baryon acoustic oscillations 



Before recombination baryons and photons are tightly coupled and gravity and pres- 
sure gradients induce sub-acoustic- Hubble-radiu s oscillations in the baryon-photon fluid 
(jSunyaev fc Zerdovichlll970l ; iPeebles fc Yulll970l ). These transmute into the acoustic peaks 
observed now in the CMB anisotropy angular power spectrum, which provide very useful in- 
formation on various cosmological parameters. The baryonic matter gravitationally interacts 
with the dark matter and so the matter power spectrum should also exhibit these "baryon 
acoustic" wiggles. Because the baryonic matter is a small fraction of the total matter the 
amplitudes of the BAO wiggles are small. The BAO peak length scale is set by the sound 
horizon at decoupling, ~ 10 2 Mpc, and so detecting the BAO peak in a real space correla- 
tion function requires observationally sampling a large volume. The B AO peak in the galaxy 



correlat ion function has recently been detected b y using SPSS da ta (jEisenstein et al.l 12005 



also see lHutsill2006l ) and by using 2dF GRS data (jCole et a 



sions of the observatio nal situation see iBlake et al. 
Percival et all j2007al fbh . 



(2007) 



2005). For more recent discus- 



Padmanabhan et al.l (120071 ). and 



The sound horizon at decoupling can be computed from relatively well-measured quan- 
tities by using relatively well-established physics. Consequently it is a standard ruler and 



can be used t o trace the universe's expansion dynamics (s ee, e.g., IBlake &: Glazebrookil2003 



Linder 2003; Seo fc Eisenstein 2003 



Hu fc Haimanl 120031 . and references therein). A mea- 



surement of the BAO peak length scale at redshift z fixes a combination of the angular 
diamet er distance and Hubbl e parameter at that redshift. More precisely, what is deter- 
mined (jEisenstein et al.l 120051 ) is the distance 



D v (z) = [(1 + z) 2 d 2 A (z)z/H(z 



1/3 



(4) 



3 The galaxy cluster gas mass fraction test was proposed bv ISasakil (|1996l ) and iPenl y.997) 
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where H(z) is the Hubble parameter and the angular diameter distance 



(1 + z)d A (z) = < 



H Q a Q Rsmhlj dz' / (H(z')a R) 
H [ dz'/H(z') 



H a Rsm[ j dz' / (H(z')a R) 



open model, 
flat model, 
closed model, 



(5) 



in the notation of iPeebled (119931 . Chap. 13). Here a^R is the radius of curvature of spatial 
hypersurfaces at fixed time and H is the present value of the Hubble parameter. Dy(z) 
depends on the cosmological parameters of the model, including those which describe dark 
energy, so we can constrain these parameters by comparing the predicted Dy(z) to the 
measurements. 



3. Computation 

In this paper we study the </>CDM model and compare our results with those derived 
in the standard ACDM model and the XCDM parametrization. In all three cases two 
parameters completely describe the background dynamics. For ACDM this pair is the non- 
relativistic matter density parameter Q m and the cosmological constant density parameter 
Qa, both defined relative to the critical energy density today. In the ACDM model we study, 
the spatial curvature density parameter £\ = 1 — fl m — £l\ need not vanish. For 0CDM 
the model parameters are Q m and a nonnegative constant a which characterizes the scalar 
field potential energy density V(0) oc (p~ a . In the 0CDM case we consider only spatially-flat 
models and 0CDM at a = is equivalent to ACDM with the same Q m and 1\ = 0. The 
XCDM parametrization is characterized by fl m and the negative equation of state parameter 
co> x = p/p. Spatial curvature is also taken to be zero for the XCDM case and XCDM at 
co> x = — 1 is equivalent to ACDM with the same f2 m and = 0. 

To compare theoretical predictions with observations we have to compute the angular 
diameter distance as a function of redshift for all three dark energy models. The angular 
diameter distance (equation [5]) depends on the Hubble parameter as a function of redshift. 
The Hubble parameters in the ACDM and XCDM models are given by 

H(z) = H ^n m (l + zf + fi A + (1 - On - O0(l + z)\ (6) 

and 
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H(z) = H ^n m {l + zf + (1 - a m )(l + z)3d+-x). ( 7 ) 

For the 0CDM model we consider H(z) does not have an analytical expression. Instead one 
has to solve the coupled system of differential equations, 



H(z) = H ^/n m (l + zf + (0 2 + k<f>~ a /G)/12, (8) 

} + 30 - ^<f>~ ia+1) = 0, (9) 

where A; is a constant. Figure 1 shows the residuals of comoving distance calculated in 0CDM 
and XCDM as a function of redshift. The 0CDM values are computed for Vt rmm = 0.3 and 
a = 2, XCDM predictions are computed for the same value of fi rmm and three diferent 
values of uj x : u x = —0.5 (solid line), u x = —1.0 (dashed line), and u x = —2.0 (dotted line). 
In general, the theoretical lines do not reduce to each other for all redshifts for any set of 
parameters a and uo x (other then a = which is the same as uo x = — 1). Figure 1 shows that 
for a = 2 the predictions of 0CDM and XCDM models already differ significantly. 

We examine the constraints on the two cosmological parameters for each dark energy 
model from two measurements of the BAO peak. The first is from the BA O peak measured 



at z = 0.35 in the correlation function of luminous red galaxies in the SDSS ( Eisenstein et al. 



20051 ). This measurement results in A(0.35) = 0.469 ±0.017 (one standard deviation error) 



where the dimensionless and ifo-independent function 

A(z) = Dv(z)^^- (10) 

and D v (z) is the distance measure defined in equation (jl]). The measured value of .4(0.35) 
does not depend on th e dark energy model a nd only weakly depends on the baryonic energy 



density (see Sec. 4.5 in lEisenstein et al.ll2005l ). The measurement also has a weak dependence 
on parameters like the spectral index of primordial scalar energy density perturbations (the 
assumed value is n = 0.98) and the sum of the neutrino masses, but this is not strong enough 
to have significant effect on the final result. To constrain cosmological model parameters in 
this case we perform a standard x 2 analysis. 

The second BAO peak measurement we use is from the correlation function of galaxy 
samples drawn from the SDSS and 2dFGRS at two different redshifts, z = 0.2 and z = 0.35, 
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as determined by iPercival et al. j2007ah FI This measurement gives the correlated values 
r s /D v (0.2) = 0.1980 ± 0.0058 and r s /D v (0.35) = 0.1094 ± 0.0033 (one standard devia- 
tion errors), •where r p is the comoving sound horizon at recombination, equation (8) of 
Percival et al.l (j2007al ). These two measurements are correlated, with the inverse of the 
correlation matrix given by 



V 



35059 -24031 
-24031 108300 



To compute r s we first compute the angular diameter distance to the surface of last scattering, 
g?a(1089). We then use the WMAP measureme nt of the apparent acoustic horizon angle in 
the CMB anisotropy data (jSpergel et al.l 120071 ) to determine the sound horizon r s = [(1 + 
z)cIa(z)]\z=w89 x 0.0104 (where we ignore the WMAP measurement uncertainty and assume 
that r s is known perfectly). The use of the WMAP prior on the apparent acoustic horizon 
angle results in very tight c o nstrain ts on the spatial curvature. When this measurement is 
not used, the IPercival et al.l (j2007al ) measurements al one can not tightly co nstrain the dark 
energy parameters (see the shaded areas in Fig. 12 of IPercival et al.ll2007al ). 



To constrain cosmological parameters in this case we follow IPercival et al.l (j2007al ) and 
first compute 

r s /D v (0.2,tt m , a) -0.1980 
r s /£V(0.35,fi m , a) -0.1094 



a) 



11) 



where for definiteness we consider the 0CDM model. We then compute the \ 2 function 



X 2 (tt m ,a) = X- 1 V- l X. 



and the likelihood function 



(12) 



L(Q m , a) = exp(-x 2 (^ m , ol)/2). 



(13) 



For both the lEisenstein et al.l (120051 ) and the IPercival et al.l (j2007al ) measurements and 



for all the models we consider, ACDM, 0CDM, and XCDM, x 2 is a function of two pa- 



rameters, either (£] m , Oa)> (^r. 



a 



or (O r 



To define lcr, 2<r, and 3a confidence level 



This analysis includes the SDSS luminous red galaxies, so the IPercival et all (|2007al ) and lEisenstein et al 



1 2005 ) BAO peak me asurements are not statistically independent. There are a number of ways to use the 
Percival et alj (|2007aT ) BAO peak measurements to constrain cosmological parameters. Here we use their 
Sic/ Dy method to compute constraints. 
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contours in these two-dimensional parameter spaces, we pick sets of points with \ 2 values 
larger than the minimum \ 2 value by 2.3, 6.17, and 11.8, respectively. 



Figures 2, 3, and 4 sh ow constraints on ACD M, XCDM, and 0CDM from the lEisenstein et al. 



( 120051 . dashed lines) and iPercival et al.l (j2007al . solid lines) data. Parts of the contours are 
not smooth because of computational noise. The BAO peak contours in these figures show 
that the measurement essentially constrains only one free parameter. When BAO peak mea- 
surements at other redshifts become available in the future, BAO data should then constrain 
both cosmological parameters. 

For now, to break this degeneracy and constrain both free parameters we use these BAO 
results together with con straints from galaxy cluster gas mass fraction ve rsus redshift data 
( ISamushia fc Ratrall2008l ). The new galaxy cluster gas mass fraction data (lAUen et al.ll2008l ) 
gives the ratio of X-ray emitting hot baryonic gas mass to total gravitational mass for 42 
hot, dynamically relaxed galaxy clusters in a redshift range from z = 0.05 to z = 1.1. Since 
the gas mass fraction of these relaxed clusters is expected to be independent of redshift these 
measurements can be used to constrain cosmological model parameters. 

In a given cosmological model, the predicted cluster gas mass fraction also depends on 
the value of the Hubble constant and the density of baryonic matter. We treat these as "nui- 
sance" parameters, assume prior probability distribution functions for them, and marginalize 
over them to derive the prob ability distribution function for the pairs of cosmological pa- 



rameters of interest (see, e.g.. lGanga et al.lll997l ). Since there still is some uncertainty in the 



values of these parameters we use two sets of Gaussian priors in our c omputations. One i s 
the set h = 0.73 ± 0.03 and VL h h 2 = 0.0223 ± . 0008 from WMA P data (jSpergel et al.ll2007h . 



the second is h = 0.68 ±0.04 ( Gott et al.ll200ll ; IChen et a 



20031) and Q h h 2 = 0.0205±0.0018 



(IFields fc Sarkai 120061). all one standard deviation errorsjj Confidence level contours derived 
from lAllen et al.l (120081 ) cluster gas mass fraction data are shown in Figs. 2, 3, and 4 as two 
sets of dotted lines corresponding to the two sets of priors for the Hubble constant and 
baryonic matter mass density. 

Since the gas mass fraction and BAO peak measurements are statistically independent 
we define the joint x 2 function by adding together the individual x 2 s- The resulting joint 
constraints are shown in Figs. 5, 6, and 7. 



5 The priors from the WMAP data that we use are derived for the best fit ACDM model. For other models 
and different parameter values the estimates on h and fib will be slightly different. Since the joint likelihood 
function does not strongly depend on h and f2b priors (for reasonable ranges of the priors) this does not lead 
to big differe nces in the final result . If one wants to be more rigor ous, slightly broader p riors from the HST 
Key Project (jFreedman et al.ll2000h and big bang nucleosynthesis (jKirkman et al.l 120031 ) can be used. 
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Results and discussion 



The lEisenstein et all ( 120051 ) BAO peak measurement has been used in co nj unction with 
other data to place constraints on various cosmological models (see, e.g.. lAlam &: Sahni 



2006 



Shafieloc 



ol 120071). 



The more recent 



Nesseris fc Perivolaropul os 2007 1: iMovahed et aljgOpJ ; [Zhang k Wull2007l : IWrightll2007l : 



purpose ( Ishida et al. 2008 ; Lazkoz et al. 20081 ). 



Perciyal et al.1 (j2007al ) data has also been used for this 



Constraints from BAO peak measurements and galaxy cluster gas mass fraction data are 
shown in Figs. 2, 3, and 4. The solid line contours in Figs. 2 and 3 show the constraints on 
ACDM and XCDM derived from the iPercival et al.1 (]2007al ) BAO data and are comparable 
to those shown with dashed lines in Fig. 12 in their paper . The d ashed contours in Fig. 3 are 
comparable to those shown in Fig. 11 in lEisenstein et al.l (120051 ). lEisenstein et al. (boosh do 
not show contours for ACDM (see our Fig. 2) and the B AO contours we show i n Fig. 4 have 
not previously been presented. Figure 2 shows that the IPercival et al.l (j2007al ) constraints, 
which make use of the WMAP measurement of the apparent acoustic horizon angle, constrain 
the sum of parameters fl\ and Q m to be very close to one (Ok = 1 — fi m — ~ 0) and favor 
a close to spatially flat model if dark energy is time independent. The spatial curvature is 
constrained so well mainly because we use the WMAP measurement of the apparent acoustic 
peak angle. BAO measurements by themselves ca n not effectively cons train dark energy 
parameters very well (see shaded areas in Fig. 12 of IPercival et al.ll2007al ). In spatially-flat 
models BAO peak measurements put tight constraints on the Q m parameter; they do not 
well constrain the "orthogonal" cosmological parameter Q\ and in particular they allow dark 
energy to vary in time (see Figs. 3 and 4). 

The BAO constraints are significantly tigh ter than the Hubble parameter versus red- 
shift dat a ones (see, e.g., Samushia et al.l 120071 ) and the strong gravitational lensing ones 
(see, e.g.. lChae et al.ll2004l ). They are, in general, about as constraining as the SNIa results 



and constrain roug hly the same linear combination of cosmol ogical parameters (see , e.g., 
Wilson et al.ll2006l ). In Figs. 2 and 3 the best fit values from IPercival et al.l (j2007al ) mea- 
surement and WMAP prior are more then 3a away from the best fit of the cluster gas mass 
fraction constraints. This is most probably due to unknown systematic errors in one or both 
of the measurements or an effect of poor statistics and should change when more and better 
data are available. 

The joint BAO peak and cluster gas mass fraction constraints are shown in Figs. 5, 6, 
and 7. They are fairly restrictive and favor a spatially-flat ACDM model with Q m ~ 0.25 
and a < 0.5 on la confidence level. Since the predictions of 0CDM for a very small value 
of the a parameter are very close to the predictions of the spatially-flat ACDM model, cur- 
rent observational tests are unable to discriminate between a time-independent cosmological 
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constant and a slowly varying scalar field with a of order 1. All three models considered 
here give about the same \ 2 ~ 52 for 41 degrees of freedom and there is no reason to favor 
one model over another based on Bayesian statistics. The constraints from the joint analy- 
sis on all three dark energy models are comparabl e to the constrain ts derived from a joint 
analysis (IWilson et al.l 120061) of earlie r SNIa data (IRiess et al.l 120041 ) and earlier cluster gas 
mass fraction data (lAUen et al.l 12004 ). Constraints on a derived from the joint analysis are 
stronger than the results quoted in previously published papers (see Table 1) The joint BAO 
peak and gas mass fraction data constraints on ACDM and XCD M derived here are a little 
weaker than t hose derived from BAO peak and more recent SNIa (lAstier et al.l 120061 ) data, 
see Figs. 13 of iPercival et al.l (j2007al ). In the joint analysis done here the uncertainties on h 
and Qbh 2 play a less significant role than they do in the cluster gas mass fraction analysis, 
i.e., the contours for the two prior sets are closer to each other in Figs. 5, 6, and 7, than in 
Figs. 2, 3, and 4. The contours in Figs. 5 and 6 are i n agreement with tighter joint results 
from other data sets considered by IWang et al.l (120071 ). 

From Fig. 1 it is clear that for large values of a 0CDM and XCDM models predict 
different cosmological evolution. For small values of a, however, if parameters are cho- 
sen appropriately, different dark energy models will at low redshifts predict very similar 
background evolution. Because of that, low redshift distance measurements have to be com- 
plemented with high redshift CMB and large scale structure measurements to discriminate 
between dark energy models. Better quality BAO peak data at a number of redshifts and 
more gas mass fraction measurements, along with tighter priors on nuisance parameters like 
the Hubble parameter and the density of baryonic matter, will allow for tighter constraints 
on dark energy parameters and could soon either detect a time dependence in dark energy 
or constrain it to a very small value. 



We acknowledge helpful discussions with M. Vogeley and support from DOE grant DE- 
FG03-99EP41093 and INTAS grant 061000017-9258. 
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Table 1. Constraints on the 0CDM a parameter from different observational tests. 



3(7 constraints on a 


Observational test(s) used 


reference 


a not well constrained 


SNIa 


Podariu & Ratra f2000) 


a not well constrained 


Radio galaxies 


Podariu et al. (2003) 


a not well constrained 


Gravitational lensing 


Ohae et al. (2004) 


a < 6.5 


Galaxy clusters 


Chen & Ratra (2004) 


a < 5 


Radio galaxies & SNIa "Gold" data set 


Wilson et al. (2006) 


a < 5 


Radio galaxies & SNIa 


Dalv et al. (2007) 


a < 4.5 


Galaxy clusters 


Samushia & Ratra (200£ 


a < 3.5 


Galaxy clusters & BAO 


This work 
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Fig. 1. — Residuals of comoving distance calculated in the 0CDM model with fi m = 0.3 and 
a = 2 and predictions of XCDM parameterization with u x = —0.5 (solid line), uj x = —1.0 
(dashed line), and u x = —2.0 (dotted line). Comoving distance is normalized to a c/H . 




Fig. 2. — 1, 2, and 3<r c onfidence level contou rs for the ACDM model. Solid lines are the 
constraints derived from iPercival et al.l (j2007al ) BAO measurements together with WMAP 
measurement of acoustic horizon an gle (the x for th e best fit value near Q m is obscured) 
and dashed lines are those from the lEisenstein et al.l (120051 ) BAO measurement. The two 
sets of dotted lines are t he constraints derived from galaxy cluster gas mass fraction data 
( ISamushia fc Ratrall2008l ) (o for the best fit value); the thick dotted lines are derived using 
the WMAP priors for h and Q^h 2 while the thin dotted lines are for the alternate priors (see 
text). The dot-dashed line corresponds to spatially-flat ACDM models. 
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Fig. 3. — 1, 2, and 3a confidence level contours for the XCDM parametrization. Solid lines 
are the constraints derived from iPercival et al.l (J2007aj) BAO measurements together with 
WMAP measurem ent of acoustic horizon angle (x for the best fit value) and dashed lines 
are those from the lEisenstein et al.l (120051 ) BAO measurement. The two sets of dotted lines 



are t he constraints derived from galaxy cluster gas mass fraction data (jSamushia fc Ratra 



20081 ) (o for the best fit value); the thick dotted lines are derived using WMAP priors for h 
and Q^h 2 while the thin dotted lines are for the alternate priors (see text). The horizontal 
dot-dashed line corresponds to spatially-flat ACDM models. 
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Fig. 4. — 1, 2, and 3<r c onfidence level contou rs for the 0CDM model. Solid lines are the 
constraints derived from iPercival et al.l (J2007aj) BAO measurements together with WMAP 
measure ment of acoustic horizon angle (o for the best fit value) and dashed lines are those 
from the lEisenstein et al.l (120051 ) BAO measurement. The two sets of dotted lines are the 
constraints derived from galaxy cluster gas mass fraction data (jSamushia &: Ratrall2008l ) (x 
for the best fit value); the thick dotted lines are derived using WMAP priors for h and Q^h 2 
while the thin dotted lines are for the alternate priors (see text). The a = axis corresponds 
to spatially-flat ACDM models. 
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Fig. 5. — 1, 2, and 3a con fidence level contours f or the ACDM model. Solid lines are the joint 
constraints derived from iPercival et al.l (J2007aj) BAO measurements together with WMAP 
measurement of acoustic horizon angle and galaxy cluster gas mass fraction data (o shows 
the best fit value with an acceptable x 2 ~ 5 7 for 42 degrees of free dom); dashed lines are 
the corresponding joint constraints using the lEisenstein et al.l (120051 ) BAO measurement (x 
shows the best fit value with an acceptable x 2 — 52 for 41 degrees of freedom). Thick lines are 
derived using the WMAP priors for h and fl^h 2 while thin lines are for the alternate priors. 
The joint best fit values for the two prior sets overlap. The dot-dashed line corresponds 
to spatially-flat ACDM models. The Q m and Q\ ranges shown here are smaller than those 
shown in Fig. 2. 




Fig. 6. — 1, 2, and 3a confidence level contours for the XCD M parametrization. Solid lines 
are the joint constraints derived from iPercival et al.l (j2007al ) BAO measurements together 
with WMAP measurement of acoustic horizon angle and galaxy cluster gas mass fraction 
data (o shows the best fit value with an acceptable \ 2 — 56. 5 for 42 degrees of free dom); 
dashed lines are the corresponding joint constraints using the lEisenstein et al.l ( 120051 ) BAO 
measurement (x shows the best fit value with an acceptable \ 2 — 52 for 41 degrees of 
freedom). Thick lines are derived using the WMAP priors for h and Q^h 2 while thin lines 
are for the alternate priors. The joint best fit values for the two prior sets overlap. The 
horizontal dot-dashed line corresponds to spatially-flat ACDM models. The Q m and u x 
ranges shown here are smaller than those shown in Fig. 3. 
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Fig. 7. — 1, 2, and 3a con fidence level contours for the 0CDM model. Solid lines are the joint 
constraints derived from iPercival et al.l (J2007aj) BAO measurements together with WMAP 
measurement of acoustic horizon angle and galaxy cluster gas mass fraction data (o shows 
the best fit value with an acceptable x 2 — 58 for 42 degrees of free dom); dashed lines are 
the corresponding joint constraints using the lEisenstein et al.l (120051 ) BAO measurement (x 
shows the best fit value with an acceptable x 2 — 52 for 41 degrees of freedom). Thick lines 
are derived using the WMAP priors for h and Q^h 2 while thin lines are for the alternate 
priors. The joint best fit values for the two prior sets overlap. The a = axis corresponds 
to spatially-flat ACDM models. The Q m and a ranges shown here are smaller than those 
shown in Fig. 4. 



